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1. 


INTRODUCTION 


This  contract  involves  research  in  support  of  three  Phillips  Laboratory  flight 
experiments,  PASP  Plus,  SPREE,  and  CHAWS. 

1.1  PASP  PLUS 


The  objective  of  the  PASP  Plus  analysis  portion  of  this  contract  is  to  examine  the 
interactions  of  the  PASP  Plus  spacecraft  and  test  solar  cells  with  the  low-energy  plasma 
environment.  Using  calculational  and  experimental  data,  we  developed  validated 
models  of  the  parasitic  current  collection  by  the  test  arrays. 

The  primary  contributors  to  this  portion  of  the  contract  were  Dr.  Victoria  A.  Davis  and 
Ms.  Barbara  Gardner. 

During  the  first  year,  we  gathered  information  on  the  PASP  Plus  arrays  including  a  site 
visit  to  Orbital  Sciences  Corporation. 

We  identified  the  major  physical  interactions  that  control  the  collection  of  parasitic 
current  from  the  space  plasma  environment  during  flight.  The  equilibrium  potential  of 
the  insulating  surfaces  exposed  to  the  plasma  depends  on  geometric  shielding  from  the 
plasma,  barrier  shielding  from  the  plasma,  and  the  orientation  of  the  surface  with 
respect  to  the  ram.  Barrier  formation  may  play  a  role  at  the  lower  biases,  particularly  for 
module  3.  Barrier  formation  occurs  at  higher  biases  under  laboratory  conditions  than 
under  space  conditions.  Snapover  plays  a  major  role  in  determining  the  current 
collected. 

We  examined  the  PASP  Plus  floating  potential.  In  order  to  determine  the  floating 
potential,  we  need  to  characterize  the  electron  collection  of  the  biased  test  solar  array, 
the  electron  collection  of  the  power  solar  array,  and  the  ion  collection  of  the  spacecraft 
body.  The  ion  collection  of  the  spacecraft  body  was  characterized  using  the 
NASCAP/LEO  computer  code.  Chapter  2  of  PL-TR-94-2212,  ADA295230,  Scientific 
Report  No.  1,  describes  this  work. 
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We  modified  the  EPS  AT  electron  density  model  so  that  it  can  model  the  1000  km  to 
2000  km  regime  that  IRI-86  does  not  properly  model.  Chapter  3  of  PL-TR-94-2212, 
ADA295230,  Scientific  Report  No.  1  describes  this  extension. 

During  the  first  and  second  years,  using  the  Gilbert  code,  we  examined  in  detail  the 
current  collection  properties  of  a  single  interconnect  or  solar  cell  gap.  We  developed  a 
formula  for  all  of  the  solar  array  geometries  that  were  flown  that  gives  the  current 
collected  as  a  function  of  the  plasma  and  electrical  parameters.  The  calculational 
technique  has  been  automated.  We  adjusted  our  computational  algorithms  to  include 
some  additional  effects  that  examination  of  the  experimental  results  suggest  are 
important.  A  discussion  of  this  work  is  in  Chapter  2  of  PL-TR-95-2095,  AD  A301487, 
Scientific  Report  No.  2. 

We  modeled  the  experiment  panels  of  PASP  Plus  to  look  for  any  interference  of  sheaths 
and  the  range  of  plasma  and  electrical  parameters  at  which  interference  might  affect  the 
experiment  results.  A  discussion  of  this  work  is  in  Chapter  3  of  PL-TR-95-2095, 
ADA301487,  Scientific  Report  No.  2. 

We  modified  the  solar  array  module  of  the  EPSAT  computer  code  for  use  by  this 
project.  The  current  to  a  solar  array  is  now  computed  by  adding  together  the  current  to 
all  the  interconnects,  gaps,  and  edges  of  the  solar  cells.  Information  on  the  environment 
is  used  by  the  formula  the  computes  the  current  to  each  interconnect,  gap,  or  edge. 
Updated  documentation  is  included  in  Chapter  12  of  PL-TR-95-2095,  ADA301487, 
Scientific  Report  No.  2. 

During  the  second  and  third  years,  we  examined  flight  data.  We  confined  ourselves 
primarily  to  data  collected  with  the  emitter  off  (so  that  density  measurements  are 
available)  and  with  the  solar  arrays  facing  the  ram.  We  compared  the  Langmuir  probe 
measurements  with  models.  We  examined  current  collection  as  a  function  of  bias  and 
plasma  density,  variation  within  a  measurement,  and  APEX  current  collection 
characteristics  and  floating  potential.  We  compared  the  measurements  with 
calculations,  explored  some  possible  reasons  for  the  scatter  in  the  data,  and  examined 
current  collection  under  negative  bias,  and  collection  by  the  APSA  array.  These 
investigations  are  discussed  in  Chapters  4  through  11  of  PL-TR-95-2095,  ADA301487, 
Scientific  Report  2  for  this  contract. 
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We  wrote  a  stand  alone  computer  program  that  implements  the  parasitic  current 
collection  model  implemented  in  EPSAT.  Documentation  is  in  Chapter  2. 

We  presented  a  paper  "Parasitic  Current  Collection  in  LEO"  at  the  SPRAT  conference  at 
NASA/LeRC  in  Cleveland  OH.  (Davis  and  Gardner,  Proceedings  of  the  XII  Space 
Photovoltaic  Research  and  Technology  Conference,  NASA  CP  3278,  p.  227) 

We  presented  a  paper  "Parasitic  Current  Collection  by  Solar  Arrays  in  LEO"  at  the 
Aerospace  Sciences  Meeting  and  Exhibit  in  Reno,  NV.  (Davis  and  Gardner,  A1 AA 
95-0594) 

We  prepared  a  paper  on  Modeling  of  Parasitic  Current  Collection  by  Solar  Arrays  in 
Low  Earth  Orbit.  It  will  be  completed  and  submitted  under  a  related  PL  contract. 


1.2  SPREE 


The  objective  of  the  SPREE  analysis  portion  of  this  contract  is  to  use  comparisons  with 
flight  data  to  examine,  improve,  and  validate  sheath  and  electron  collection  models  and 
to  examine  some  of  the  interactions  that  determine  the  spectrum  of  the  incident  ions. 

The  primary  contributors  to  this  portion  of  the  contract  were  Dr.  Gary  A.  Jongeward, 
Dr.  Ira  Katz,  and  Dr.  Victoria  A.  Davis. 

During  the  first  year,  we  developed  a  module  for  EPSAT  that  examines  effluent 
ionization  for  shuttle  thrusters.  We  defined  860  lb  thrusters  and  computed  the 
ionization  due  to  the  interaction  of  the  electron  beam  and  the  thruster  effluent.  We 
added  a  module  that  computes  the  density  of  ions  created  by  effluent  ionization  as  a 
function  of  position.  We  developed  a  detector  module  that  computes  the  flux  of  ions 
created  by  effluent  ionization  to  the  detector.  Chapter  4  of  PL-TR-94-2212,  AD  A295230, 
Scientific  Report  No.  1  describes  this  work. 

During  the  second  year,  we  used  the  EPSAT  model  to  compare  the  computed  spectrum 
with  SPREE  measurements  for  a  period  during  which  the  FPEG  beam  was 
perpendicular  to  the  Earth's  magnetic  field.  This  is  the  event  described  in  Hardy,  et  al. 
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(Observations  of  Electron  Beam  Propagation  Perpendicular  to  the  Earth's  Magnetic 
Field  during  the  TSS 1  Mission).  A  discussion  of  the  calculation  is  in  Chapter  13  of 
PL-TR-95-2095,  Scientific  Report  No.  2  and  documentation  for  the  model  is  in  Chapter 
14  of  PL-TR-95-2095,  ADA301487,  Scientific  Report  No.  2. 


During  the  third  year,  we  worked  on  a  comparison  of  the  shuttle  floating  potential  as 
meastured  by  SPREE  with  the  Floating  potential  determined  from  the  magnetic  field, 
velocity,  and  tether  current.  This  work  is  described  in  Chapter  3. 

We  estimated  the  LC  time  constant  for  the  shuttle  charging  process  to  be  15  Hz.  This 
does  not  provide  a  physical  explaination  for  the  1  Hz  modulation  of  the  ram  flux  seen  in 
daylight  by  SPREE  under  some  conditions. 

We  assisted  Phillips  Laboratory  staff  with  a  paper  on  thruster  induced  charging. 

The  paper  "Observations  of  Ionosphere  Heating  in  the  TSS-1  Subsatellite  Presheath" 
was  published  in  the  May  1, 1994  issue  of  Jouvnul  of  Geophysical  Rcseavch.  (Katz,  et  al, 
]GR,  99,  p.  8961). 

We  presented  a  paper  "Observations  of  Ionosphere  Heating  in  the  TSS-1  Subsatellite 
Presheath"  at  the  4th  International  Conference  on  Tethers  in  Space  in  Washington,  DC. 


1.3  CHAWS 


The  objective  of  the  CHAWS  analysis  portion  of  this  contract  is  to  develop  and  validate 
models  of  the  spacecraft  wake,  including  the  currents  to  the  various  detectors  of  the 
CHAWS  wake  side  sensor  on  the  Wake  Shield  Facility. 

The  primary  contributors  to  this  portion  of  the  contract  were  Dr.  Victoria  A.  Davis  and 
Dr.  Myron  M.  Mandell. 

During  the  first  year,  we  used  our  two-dimensional  computational  tools  to  compute  the 
current  to  a  probe  behind  a  disk  in  the  geometry  used  in  the  Hanscom  laboratory  tank. 
Chapter  5  of  PL-TR-94-2212,  AD  A295230,  Scientific  Report  No.  1  describes  this  work. 
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We  used  our  two-dimensional  computational  tools  to  estimate  the  relative  importance 
of  the  thermal  and  electrostatic  effects  for  the  flight  experiment.  Chapter  6  of 
PL-TR-94-2212,  ADA295230,  Scientific  Report  No.  1  describes  this  work. 

We  used  DynaPAC  to  compute  the  current  collected  by  the  CHAWS  wake  side 
detectors  for  a  range  of  flight  conditions.  We  developed  a  scheme  for  interpolating  to 
additional  conditions.  We  examined  the  effect  on  current  collection  of  the  high  voltage 
screen  and  grounded  current  collection  plate  design.  Chapter  7  of  PL-TR-94-2212, 
ADA295230,  Scientific  Report  No.  1  describes  this  work. 

We  developed  a  screen  dipole  model  that  reproduces  some  of  the  qualitative  features  of 
the  three-dimensional  calculations.  Chapter  8  of  PL-TR-94-2212,  ADA295230,  Scientific 
Report  No.  1  describes  this  model. 

We  developed  two  new  modules  for  the  EPSAT  computer  code.  The  Wake  Studies 
module  explores  the  density  and  current  to  a  detector  in  the  ambient  plasma  or  behind 
a  disk  in  the  neutral  approximation.  This  module  was  revised  through  the  contract.  The 
final  documentation  is  in  Chapter  4.  The  new  CHAWS  Current  Module  gives  the 
current  to  the  wake  side  probe  as  computed  by  the  D5maPAC  and  POLAR  computer 
codes.  Chapter  10  of  PL-TR-94-2212,  Scientific  Report  No.  1  describes  this  module. 

The  geomagnetic  coordinates  were  added  to  the  Magnetic  Field  at  a  Point  screen  of 
EPSAT.  Chapter  11  of  PL-TR-94-2212,  ADA295230,  Scientific  Report  No.  1  describes  the 
implementation. 

The  EPSAT  floating  potential  model  was  modified  to  include  the  effect  of  the  variation 
in  the  v  x  B  potential  across  an  object.  Chapter  12  of  PL-TR-94-2212,  ADA295230, 
Scientific  Report  No.  1  describes  this  implementation.  Two  new  subscreens  to  the 
floating  potential  screen  were  added.  One  shows  the  minimum  and  maximum 
potentials  at  the  comers  of  each  object.  The  other  allows  the  user  to  display  the 
"surface"  potential  on  the  entire  system. 

We  used  the  measurements  of  the  secondary  emission  yield  of  stainless  steel  made  at 
Hanscom  as  a  basis  for  a  formula  for  the  secondary  emission  as  a  function  of  incident 
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particle  energy  and  species.  Chapter  13  of  PL-TR-94-2212,  ADA295230,  Scientific  Report 
No.  1  describes  this  work. 

During  the  second  year,  we  examined  data  collected  during  the  flight.  We  wrote  some 
simple  computer  programs  to  assist  us  in  the  data  analysis. 

The  EPSAT  model  of  WSF  and  the  shuttle  was  modified  to  obtain  better  values  of  the 
floating  potential.  Hidden  surfaces  of  the  engine  bells  that  cannot  collect  current  were 
defined  to  be  insulating  so  that  they  do  not  collect  ram  current.  The  WSF  carrier  was 
added  to  the  EPSAT  model. 

We  did  rough  calculations  of  the  fraction  of  phase  space  of  the  ram  ions  incident  on 
WSF  that  reach  the  center  of  the  wake  side  of  the  disk. 

We  presented  a  paper  "CHAWS  Wake  Side  Current  Measurements:  Comparison  with 
Preflight  Predictions"  at  the  33rd  Aerospace  Sciences  Meeting  in  Reno,  NV.  (Mandell, 
effl/,AIAA  95-0489). 
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2. 


DOCUMENTATION  FOR  SOLAR  PROGRAM 


Parasitic  current  collection  is  determined  by  die  solar  array,  the  plasma 
environment,  and  the  applied  bias.  The  user  can  input  values  that  describe  the 
array  and  the  plasma  environment.  The  program  solar  generates  an  iv-curve. 

2.1  Cell  Properties 

Acj)  =  Potential  generated  by  a  single  solar  cell. 

Vsnap  “  Snapover  Potential  (First  Crossover  Energy) 
mi  =  Ion  Current  Slope 
Tcell  Cell  Thickness 
Dgap  =  Gap  Width 

Teg  =  Coverglass  and  Adhesive  Thickness 
Tadh  “  Adhesive  Thickness 
To  =  Coverglass  Overhang 

Figure  1  illustrates  the  solar  cell  gap  region  which  is  modeled  using  the  above 
parameters. 
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Figure  1.  Solar  cell  gap  region. 

2.2  Array  Properties 

Li  =  Total  Interior  Length  without  Interconnect 
Lie  =  Total  Interconnect  Length 
Le  =  Total  Edge  Length  without  Metal 
Lm  =  Total  Edge  Length  with  Metal 

N  =  Number  of  Cells  per  Series 
Ns  =  Number  of  Series 

The  array  geometry  is  shown  in  Figure  2. 
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Figure  2.  Solar  Array  with  12  cells  in  series. 

The  relationship  between  the  cell  properties  and  the  array  properties  for  this 
array  is: 

Lie  =  Interconnect  length*number  between  cells 
Lm  =  Interconnect  Iength*ntunber  along  edges 
Le  =2NwWc  +  2NiLc-Lm 
Li  =  (Nw-1)NlLc+(Nl-1)NwWc  -  Lie 
2.3  Environment  Parameters 
rio  =  Plasma  Density 
0  =  Plasma  Temperature 
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2.4  Current  Collected 


The  current  collected  from  the  plasma  by  a  solar  cell  is  collected  on  the  sides  of 
the  solar  cell,  and  by  the  interconnects,  if  any.  It  is  calculated  as  the  product  of 
the  solar  cell  effective  collecting  area  and  the  incident  current  density  of  the 
attracted  species.  The  attracted  species  depends  on  whether  the  cell  potential 
with  respect  to  the  plasma  is  positive  or  negative. 

The  effective  area  is  a  function  of  the  electrical  and  geometric  properties  of  the 
cell  and  environmental  parameters  and  is  calculated  as  described  below. 

2.5  Single  Cell  Electron  Current  Collection 

For  electrons  the  effective  collecting  area  is  calculated  as  the  sum  of  four  terms: 

-‘■cou(^c)  W‘(l|.,)  +  Le  W»(<,,)  +  L^  Wf  (<.,)  , 


Where, 

Acoll  =  Solar  cell  effective  collecting  area 

Wg^  =  Sheath  width  for  interconnects 

Wg  =  Sheath  width  for  gaps  without  an  interconnect 

Wg  =  Sheath  width  for  edges 

(j)c  =  Solar  cell  potential  with  respect  to  the  plasma 

The  area  is  computed  from  several  two-dimensional  segments.  It  is  the  distance 
along  the  solar  cell  edge  times  the  length  of  the  two-dimensional  sheath.  The 
technique  to  calculate  the  sheath  width  is  described  in  the  section  "Solar  Cell 
Sheath  Width  Model." 
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Figure  3.  Solar  cell  gap  with  sheath. 

2.6  Single  Cell  Ion  Ciurent  Collection 

Ion  collection,  which  occurs  when  the  cell  potential  is  negative  with  respect  to 
the  plasma,  is  given  as  a  linear  function  of  the  cell  potential.  The  effective 
collecting  area  is  as  follows: 

Acoll(<l)c)  =  mi<l>c  (2) 

2.7  Solar  Cell  Sheath  Width  Model 


The  two-dimensional  electron  collecting  sheath  width  for  a  solar  cell  is  modeled 
using  the  following  formula: 


S.  W.  =  l.l  A(geom.)B(ii) 


(0.1,  Ilf 


max  (0. 1,  T|) 


0.44- 


7  1.13 

(j)  7^-'' 


where  the  function  B  is 


no  interconnect  geometry  and  r)  >  0.6 
otherwise 
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T]  =  min(l,lst  crossover  potential / cell  potential)  =  min(l,V snap/  <t>c) 

(|)  =  Cell  potential/  plasma  temperature  =  6 

X.  =  Debye  length  /  cell  thickness  =  ^So0/e  n,,  /Tceii 

This  formula  should  not  be  applied  to  scaled  potential  values  greater  than 
500/0.05. 

2.8  Solar  Array  Current 

The  solar  array  is  assumed  to  be  negatively  grounded. 

To  get  the  total  current  collected  by  the  solar  array,  we  sum  the  currents  to  all  the 
solar  cells  as  follows: 


N 

Isa  (<t>a)  =  Ns  X ^coll  (^a  +TlA(^)  (  (t>a  +nA(t)) 

n=l 

Where, 


if  (j)  >  0 

=  attracted  current 

if  ([)<  0 

(|)a=  Array  Potential 

2.9  Keyword  Input 


1th  =  eno^ 


eO 


27C  m. 


(4) 


(5) 


This  program  uses  keyword  input.  The  input  stream  consists  of  a  list  of 
keywords  and  values  as  given  in  Table  1.  The  END  keyword  signals  that  a 
current-voltage  curve  is  to  be  generated  with  the  present  values.  Any  number 
(up  to  1000)  of  sets  of  input  parameters  separated  by  the  END  keyword  may  be 
included.  Execution  stops  when  an  end-of-file  is  sensed. 
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Table  1.  Keywords  for  Input 


Keyword 

Symbol  Used 
Above 

Units 

BIAS 

Vi,V2,Npot 

V 

POTENTIAL 

A(t) 

V 

SNAP 

Vsnap 

V 

ISLP 

mi 

m2V-l 

CELL 

TceU 

m 

GAP 

Dgap 

m 

CVGL 

Teg 

m 

GLLOE 

Tadh 

m 

OVHG 

To 

m 

LNI 

Li 

m 

LI 

Lie 

m 

LE 

Le 

m 

LEM 

Lm 

m 

NUMBER 

N 

DENSITY 

no 

TEMP 

0 

eV 

IONS 

A  m"^ 

END 

Description 

First  potential,  last  potential,  and 
number  of  potentials  on  IV  curve 
Potential  generated  by  a  single  solar 
cell. 

Snapover  Potential  (First  Crossover 
Energy) 

Ion  Current  Slope 
Cell  Thickness 
Gap  Width 

Coverglass  and  Adhesive  Thickness 

Adhesive  Thickness 

Coverglass  Overhang 

Total  Interior  Length  without 

Interconnect 

Total  Interconnect  Length 
Total  Edge  Length  without  Metal 
Total  Edge  Length  with  Metal 
Number  of  Cells  per  Series 
Plasma  Density 
Plasma  Temperature 
Incident  Ion  Current  Density 
End  of  Data  Set 
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3. 


MEASUREMENT  OF  NEGATIVE  CHARGING  ION  DISTRIBUTIONS 
DURING  TSS-1 


Identification  of  peaks  in  ion  and  electron  energy  spectra  have  been  used  for 
decades  by  researchers  to  determine  charging  levels  on  spacecraft.  During  the 
TSS-1  mission,  a  computer  performed  this  analysis  in  real  time  and  provided 
orbiter  potentials  to  the  crevr  and  ground  controllers.  TSS-1  was  a  unique 
opportunity  to  test  the  automated  charging  level  analysis,  because,  for  the 
voltages  achieved  during  the  first  mission,  the  tethered  subsatellite  acted  as  a 
floating  probe,  and  provided  an  independent  measure  of  the  charging  potential. 
Below  we  present  examples  of  the  highly  anisotropic  ion  spectral  data  collected 
during  the  mission,  potentials  calculated  in  real  time  using  this  data,  and  show 
comparisons  with  potentials  derived  using  tether  current  and  voltage 
measurements. 

3.1  Anisotropic  Ion  Flux 

Figure  4a  shows  the  ion  flux  to  SPREE  as  a  function  of  bin  as  the  instrument 
rotates  for  day  218  at  2:36:02  GMT.  Figure  4b  shows  the  same  ion  flux  as  a 
function  of  incident  energy.  At  the  time  the  potential  decays,  (as  the  FPEG  is 
turned  off),  the  aperture  is  facing  the  ram  direction.  Note  that  the  bin  (angle) 
which  receives  the  peak  current  moves  from  5  to  0  as  the  particle  energy  shifts 
form  40  to  0  V.  The  anisotropic,  beaming  of  ions  seen  here  is  consistent  with  a 
local  potential  structure  around  the  instrument,  that  is  strong  electric  fields  near 
the  instrtunent,  and  not  a  uniformly  charged  orbiter. 
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Figure  4.  Ion  flux  to  SPREE  as  the  instrument  rotates  for  day  218  at  2:36:02  GMT. 
(a)  shows  the  ion  flux  as  a  function  of  bin.  (b)  shows  the  ion  flux  as  a 
function  of  incident  energy. 


We  did  a  series  of  calculations  to  determine  the  local  potential  structure  and 
typical  particle  trajectories.  Figure  5  shows  a  model  of  the  SPREE  instrument  and 
the  near  by  Langmuir  probe.  This  model  is  placed  in  the  nested  grid  structure 
shown  in  Figure  6.  The  coordinate  system  is  +X  toward  the  nose,  +Y  toward  the 
starboard  wing,  and  +Z  into  the  bay.  The  space  potentials  and  sheath  currents 
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Figure  6.  Grid  Structure  used  for  calculations. 


Figure  7  shows  the  computed  potentials  for  37.5  V.  Figure  8  shows  some  typical 
trajectories  in  the  potentials  shown  in  Figure  7. 


Z-AXIS 


COLOR  LEGEND 


X-AXIS 

Minimum  Potential  =  -3.75E+01  Maximum  Potential  =  9.32E-02 
6.00  <X<  17.00,  5.00  <Z<  14.00,  CUTPLANE  OFFSET  Y=  13.50 


CO 


6.  7.  8.  9.  10.  11.  12.  13.  14.  15.  16. 


X-AXIS 

Minimum  Potential  =  -3.75E+01  Maximum  Potential  =  -4.60E-10 
6.00  <X<  16.00,  9.00  <Y<  19.00,  CUTPL AN E  OFFSET  Z=  9.00 


Figure  7.  Computed  Potentials  for  a  ground  potential  of  37.5.  V. 
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X-AXIS  Z-AXIS 


X-AXIS 


Figure  8.  Typical  Trajectories  for  potentials  shown  in  Figure  7. 
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From  the  total  collected  current  we  sorted  out  the  current  that  would  be 
measured  by  SPREE.  We  used  Figure  9  to  determine  the  aperture  location  in  grid 
coordinates  to  be  (10.3, 13.95, 8.9). 


X-AXIS  X-AXIS 


9.00  <X<  11.00,  8.00  <2<  10.00.  CUTPLANE  OFFSET  Y=  14.00  9.00  <X<  11.00,  13.00  <Y<  15.00,  CUTPLANE  OFFSET Z=  9.00 

Figure  9.  Figures  used  to  determine  aperture  location  in  grid  co-ordinates. 

We  surveyed  all  the  macroparticles  hitting  the  screen  within  10°  of  the  normal  in 
the  X-Z  plane  and  within  10°  of  the  normal  to  the  X-Z  plane. 

Algebraically —Let  (x,y,z)  be  the  direction  from  the  point  of  incidence  on  the 
screen  to  (10.3, 13.95, 8.9).  Let  (u,v,w)  be  the  direction  of  travel  at  the  point  of 
incidence.  The  criteria  are  then 


The  macroparticles  were  then  binned  by  angle  from  vertical.  The  angle  is 
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Figure  10  shows  the  incident  flux  as  a  function  of  shuttle  potential  and  incident 
bin. 


Calculated  Current  Measured  by  SPREE  (A) 


Figure  10.  Computed  ion  flux  as  a  function  of  shuttle  potential  and  incident  bin. 


The  similarity  between  Figures  4  and  10  indicates  that  the  local  potentials  due  to 
the  shape  of  SPREE  create  the  dependence  of  incident  angle  on  shuttle  potential. 


3.2  SPREE  Measurements  During  DEPl  Operations 

During  DEPl  operations,  a  sequence  of  resistors  was  inserted  into  the  tether 
circuit  and  the  FPEG  electron  gun  was  cycled.  We  only  address  charging  data 
obtained  when  the  gun  was  off.  The  largest  resistance  (2.5  MD)  was  orders  of 
magnitude  larger  than  the  tether  resistance  (2200  Q) .  When  switched  into  the 
tether  circuit,  the  voltage  measured  at  the  orbiter  was  the  v  x  B*1  generated  by 
the  tether  moving  through  tire  earth's  magnetic  field,  since  the  subsatellite 
floating  potential  was  only  a  few  kT.  For  smaller  resistances,  the  maximum 
current  that  the  induced  EMF  could  drive  through  the  tether  resistance  was  often 
greater  than  the  ion  current  that  the  orbiter  could  collect.  When  this  occurred,  the 
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orbiter  charged  to  a  negative  potential .  The  steady  state  orbiter  potential  was 
such  that  the  sheath  drop  at  the  orbiter  reduced  the  voltage  across  the  tether 
enough  to  reduce  the  tether  current  to  the  available  orbiter  ion  current.  The 
voltage  drop  in  the  tether  is  product  of  the  tether  current  and  resistance,  both  of 
which  are  known.  The  sheath  drop,  and  therefore  the  orbiter  potential,  is  just  the 
difference  between  the  available  EMF  and  the  voltage  drop  in  the  tether.  The 
sheath  drop  at  the  subsatellite,  as  measured  by  RETE,  was  very  small  when 
FPEG  was  not  operating. 

In  Figure  11,  the  orbiter  potentials  for  11  DEPl  cycles  where  charging  occurred 
calculated  from  the  tether  current  and  voltage  measurements  are  compared  with 
the  orbiter  potential  as  reported  by  SPREE.  The  orbiter  charged  as  much  as  26  V 
negative,  about  50  percent  of  the  available  motion  induced  v  x  B.  The  largest 
negative  potentials  correspond  to  low  ionosphere  densities  where  the  available 
ram  ion  current  is  smallest.  When  the  available  current  is  small  the  resistive 
voltage  drop  in  the  tether  is  small.  The  magnitude  of  the  ram  ion  current  is 
corisistent  with  an  effective  ram  facing  exposed  conducting  area  of  about  23  m^ . 
This  number  agrees  with  preflight  calculations  performed  using  a  three  dimen¬ 
sional  spacecraft  charging  model  of  the  orbiter.  Those  calculations  show  tiiat 
most  of  the  collection  is  by  main  engine  bells.  Thus,  the  orbiter  potentials  as 
measured  ion  energies,  SPREE,  and  by  insertion  of  a  high  impedance  in  the 
tether,  SETS,  both  agree  with  the  basic  theory  that  the  limited  ion  sheath  current 
drives  the  orbiter  potential.  The  higher  the  available  ion  current,  the  less  the 
charging. 


22 


35.0 
30.0 
25.0 

> 

^  -  20.0 

®  15.0 
O 
Q. 

’  10.0 

5.0 
0.0 

Figure  11. 

The  SPREE  results,  with  their  associated  uncertainty  due  to  the  detector  energy 
bin  widths,  agree  quite  well  with  the  tether  results.  The  SPREE  potential  is 
shown  as  the  top  and  bottom  of  the  detector  energy  channel  which  contains  the 
most  counts.  The  SPREE  algorithm  reports  the  middle  of  the  energy  channel. 
Note  that  the  SPREE  results  were  obtained  in  real  time.  The  SETS  potentials  are 
based  on  tether  currents  and  voltages  measured  at  two  different  times,  and 
depend  on  the  particular  operational  sequences  and  accurate  subsatellite  sheath 
potentials.  These  measurements  will  not  be  available  at  high  potentials  (long 
tether  lengths)  and  magnetic  limiting  of  electron  collection  (low  plasma 
densities) ,  Or  when  the  DCORE  electron  gun  is  operating.  Under  those 
conditions,  only  the  SPREE  potentials  will  be  available. 

3.3  Shunter  Orbit 

A  synoptic  comparison  of  SPREE  measurements  over  an  entire  orbit  can  be  made 
with  potentials  calculated  using  the  measured  tether  current  and  a  model  of  the 
V  X  B  induced  EMF.  During  an  entire  orbit  from  218/10:45  to  218/12:15,  the 
tether  was  cormected  Ihrough  a  low  resistance  shunt  to  the  orbiter  chassis.  The 
orbiter  potentials  as  measured  by  SPREE  are  plotted  in  Figure  12.  The  plot  shows 
clearly  that  the  SPREE  resolution  is  limited  by  the  instrument  energy  bin  widths. 
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Shunter  Orbit 


Time  (sec)  on  day  218 


Figure  12.  Orbiter  potentials  as  measured  by  SPREE. 

The  tether  current,  as  measured  by  SETS,  is  shown  in  Figure  13.  Note  that  the 
orbiter  charges  when  the  tether  current  is  very  small,  which  occurs  when  the 
ambient  ionosphere  density  is  low.  The  reason  that  the  charging  levels  increase 
when  the  tether  voltage  decreases  can  be  seen  in  Figure  14.  Plotted  in  this  figure 
are  the  model  available  EMF,  the  tether  voltage  drop  (Ihe  product  of  the  mea¬ 
sured  tedier  current  from  SETS  and  the  tether  resistance),  and  the  difference, 
which  is  the  potential  drop  across  the  sheath  which  forms  between  die  charged 
orbiter  and  the  ionosphere.  The  tether  resistance  value  used  here  is  2200  n.  The 
charging  occurs  when  the  resistive  voltage  drop  (product  of  tether  current  and 
resistance)  is  small,  leaving  most  of  the  motion  induced  EMF  to  be  dropped 
across  the  electron  collecting  sheath  at  the  subsateUite,  and  the  ion  collecting 
sheath  at  the  orbiter.  For  these  low  potentials,  the  earth's  magnetic  field  does  not 
limit  the  electron  collection.  Most  of  the  potential  drop  occurs  across  the  orbiter 
ion  collecting  sheath. 
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Figure  13.  Tether  current  as  measured  by  SETS. 
Shunter  Orbit 


Figure  14.  Model  available  EMF  values,  tether  voltage  drop  (the  product  of  the 
measured  tether  current  from  SETS  and  the  tether  resistance),  and  the 
difference,  which  is  the  potential  drop  across  the  sheath  which  forms 
between  the  charged  orbiter  and  the  ionosphere. 
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This  calculation  of  the  orbiter  potential  is  compared  with  the  SPREE  result  in 
Figure  15.  Over  the  large  range  of  plasma  parameters  encountered  during  the 
orbit,  the  two  methods  of  calculating  the  orbiter  potential  are  in  agreement.  Once 
again,  however,  the  calculation  using  the  tether  current  wiU  not  be  available 
during  gun  operations  and  for  very  low  potentials.  With  the  DCORE  electron 
gun  operational  or  at  high  (kilovolt)  potentials,  only  the  SPREE  results  will  be 
available. 

Shunter  Orbit 


38000  39000  40000  41000  42000  43000  44000  45000 

Time  (sec)  on  day  21 8 


Figure  15.  Orbiter  potential  as  measured  by  SPREE  and  as  computed 
from  the  available  EMF  and  the  tether  current. 


3.4  Conclusions 

SPREE  measured  ram  ions  accelerated  in  a  local  sheath.  There  was  no  evidence 
of  locally  generated  ions.  A  direct  comparison  with  other  measurements 
confirms  SPREE  potential  measurements.  SPREE  will  provide  real  time 
potentials  at  high  potentials  when  the  subsatellite  is  no  longer  at  ionosphere 
potential. 
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4. 


EPSAT  WAKE  STUDIES  SCREEN 


The  EPSAT  Wake  Studies  screen  uses  the  neutral  approximation  to  compute  the 
densities  and  current  densities  (as  it  would  be  measured  by  an  RP A)  of  the 
neutral  and  plasma  species  in  front  of  or  behind  a  disk.  The  input  parameters  are 
the  disk  parameters,  the  spacecraft  velocity,  the  ambient  neutral  and  plasma 
environment,  the  point  of  interest,  the  look  direction,  and  the  incident  energy. 

The  disk  is  described  by  its  normal,  its  radius,  and  the  location  of  its  center.  The 
normal  of  either  the  ram  facing  side  or  the  wake  facing  side  can  be  given. 

By  default  the  spacecraft  velocity  is  given  by  the  EPSAT  orbit  generator. 

The  ambient  neutral  and  plasma  environments  are  shown  on  the  Neutral  Density 
at  a  Point  and  Plasma  Density  at  a  Point  screens.  The  density  of  each  neutral  and 
ion  species  and  the  neutral  and  ion  temperatures  are  used. 

The  point  of  interest  can  be  specified  either  by  specifying  an  origin  and  a 
distance  from  that  origin  in  spherical  coordinates  or  directly  in  Cartesian 
coordinates. 

The  look  direction  can  be  specified  eidier  in  terms  of  the  spherical  coordinate 
angles  or  in  Cartesian  coordinates.  This  parameter  is  only  used  for  the  current 
density  computation. 

The  incident  energy  is  given  in  electron  volts.  This  parameter  is  only  used  for  the 
current  density  computation.  If  the  energy  is  negative,  the  current  density  is 
zero. 


4.1  Densities 

On  the  ram  facing  side  of  the  disk,  the  densities  are  the  ambient  densities.  The 
point  of  interest  is  on  the  ram  facing  side  if  the  dot  product  of  the  normal  of  the 
wake  facing  side  of  the  disk  with  the  vector  from  the  center  of  the  disk  to  the 
point  of  interest  is  negative. 
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On  the  wake  facing  side  of  the  disk,  the  density  for  each  species  is  the  integral 


dQ  dE^E  g(r,  9,(1))  exp 


^  ^ram  ^V^E 


ram 


cos  9 


(9) 


J) 


The  ambient  density  of  the  species  is  pa-  The  ion  temperature  is  T.  The  vector  r  is 

the  vector  from  the  center  of  the  disk  to  the  point  of  interest.  The  energy  E  is  the 
incident  particle  energy.  The  energy  Eram  is  given  by 


E 


ram 


1 

=  — mcV 
2  ® 


2 

ram 


(10) 


where  ms  is  the  mass  of  the  species  and  Vram  is  the  spacecraft  velocity.  The 
integral  is  done  over  the  polar  coordinates  9  and  (j),  representing  die  negative  of 
the  incoming  particle's  velocity  direction,  where  the  origin  for  the  integral  is  at 
the  point  of  interest  and  9  is  measured  from  the  negative  spacecraft  velocity 
direction.  The  function  g  is  zero  if  the  direction  given  by  9  and  ([)  is  blocked  by 
the  disk  and  unity  otherwise. 

The  integral  is  solved  differently  for  the  case  where  the  disk  is  normal  to  the 
spacecraft  velocity  and  where  it  is  not.  A  different  technique  is  also  used  for  the 
case  in  which  the  point  of  interest  is  directly  behind  the  disk  and  where  it  is  to 
the  side. 

4.2  Solution  Technique  for  Normal  Orientation  and  Point  Directly  Behind 
Disk 


The  outermost  integral  over  the  azimuthal  angle,  (j),  is  performed  taking  zero  as 
the  direction  to  the  nearest  disk  edge  point.  To  get  a  concentration  of  points  near 
the  nearest  edge  point,  we  define  \|/  =  The  minimum  value  of  cos9  can  be 

shown  to  be 


t 


Umin(<t>)=-[l+(x/z)^ 


1/2 


x/r  =  -cos(t)  + 


cos^(l)+(R/ r)^-l 


1/2 


(11) 
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(The  maximum  value  of  cos0  is  1,  for  particles  overtaking  the  spacecraft  from 
behind.)  The  energy  variable  is  transformed  by  t  =  ^  where  Eq  is  an 

enerutational  convenience.  It  is  taken  to  be 


Eo 


V  +  5. 


2eT 
'  m 


In  (n.) 


(12) 


where  nt  is  the  number  of  integration  points  used  for  the  t-variable  (i.e.,  energy). 
(Note:  Eo  should  be  set  to  at  least  the  ion  temperature,  T.) 


After  appropriate  manipulation,  the  density  at  a  point  in  the  wake  of  the  disk  is 
given  by 

=  - e-^ram/T  jjdtK(i|;,t) 

V  ^  Eram  ^  (1 

K  (\|/,  t)  =  tl^^  je-YltlUminlV)  _g-Y(t)u„,4x|/)  j 


Y(t)  =  . 


-4Eo  E^Jnt 


Note  that,  to  avoid  overflows,  all  the  exponents  (including  the  one  in  the 
prefactor)  should  be  grouped  together,  so  that  the  only  exponentiations  are  the 
two  in  the  inner  loop. 


EPSAT  performs  the  integrals  numerically. 

4.3  Solution  Technique  for  Normal  Orientation  and  Point  Not  Directly 
Behind  Disk 


The  integral  is  over  those  trajectories  blocked  by  the  disk,  and  the  result 
subtracted  from  unity.  The  outer  integral  is  over  the  angle  between  the  look 
direction  and  the  vector  toward  the  ram,  which  takes  the  range 
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■<u  =  cos0<- 


To  concentrate  points  toward  the  ram  direction,  the  integration  variable  is 
changed  to  w  =  Vl-u .  The  integral  over  azimuthal  angles  contributes  a  factor 


A  (j)  (u)  =  2  arccos 


[(r^  -R^  -z^)u^  +z^j 
[2rzu(l-u’)'"]  1 


The  final  result  is 


-|3/2 

1_A=  ^  f 


Po  L^T 


r  V  min  - 

wdwA(j)(l-  (l— w^,t) 


K'(u,t)  =  t'^  ^  (-/nt)^^^ 


/  \  1-4  Eg  E__/nt 

y(‘)=J - 


4.4  Solution  Technique  for  Other  Orientations 

If  the  normal  to  the  disk  is  not  parallel  to  the  velocity,  the  problem  is  solved 
similarly,  although  the  geometry  is  more  complex. 

First,  the  disk  and  the  field  point  are  projected  onto  a  plane  normal  to  the 
velocity  (Z),  so  that  the  disk  edge  is  an  ellipse  with  its  major  axis  in  the  X 
direction. 


x^  +ay^  =R^ 
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The  point  (x,  y)  on  the  ellipse  nearest  the  field  point  projection  (xq,  yo)  (for  the 
first  quadrant)  satisfies 

=  0 

(18) 

When  (xQ,  yo)  lies  inside  the  ellipse,  the  outer  angular  integral  is  over  the  angle,  (j), 
from  the  field  point  to  the  wake  edge,  given  by 

<|)  =  arctan  — — ^  ±\|/^ 

[x-xoj  (19) 

where  the  integration  variable,  \|/,  concentrates  near  the  presumed  integrand 
peak  at  (x,  y).  There  are  two  values  of  Umin  (the  cosine  of  the  angle  between  the 

direction  from  the  field  point  to  the  disk  edge  and  the  negative  velocity 
direction)  for  each  V|/,  obtained  by  calculating  the  disk  edge  point  (xq  +  d  cos(j),  yo 
+  d  sine]),  zi).  The  values  are  determined  by  first  solving  the  ellipse  equation  d, 
then  the  circle  equation  for  zi.  Then 

Zq-Zj 

^mm 

where  the  sign  of  zi  is  opposite  the  sign  of  yo  +  d  sinq.  The  energy  integral  then 
proceeds  as  for  the  velocity-aligned  case,  except  that  the  normalization  must  take 
account  of  the  asymmetry  in  the  (j)  integral. 

When  (xo,  yo)  lies  outside  the  ellipse,  the  integral  is  over  the  region  from  which 
ions  are  excluded  and  the  normalized  result  subtracted  from  unity.  The  outer 
angular  integral  is  over  the  variable  u  (the  cosine  of  the  angle  between  the  ram 
direction  and  the  direction  from  the  field  point  to  a  point  on  the  disk),  which 
goes  from  the  angle  corresponding  to  the  near  point  of  the  disk  to  the 
diametrically  opposite  point.  For  each  u,  there  are  two  values  of  (j)  that  satisfy 
both  the  disk  edge  equation  and  the  defining  equation  for  u: 


x(a  - 1)  -  axo  +  xyo 


'r2-x2 
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(xj  +  d  cos  (j))^  +  a(yo  +  d  sincj))^  = 

(xo  +  d  cos  (|))^  +  (yo  +  d  sm(j))^  +  z^  =R^ 
z(yo  +  sin  ({))  <  0 

Z-Zq 
u  = 

The  difference  between  the  two  (j)  values  corresponds  to  the  A(|)(u)  that  appears  in 
the  velocity-aligned  case. 

4.5  Current  Densities 


The  differential  ram  ion  current  density  to  a  detector  in  the  presence  of  a  disk  is 
zero  if  the  detector  is  looking  at  the  disk.  Otherwise  it  is  given  by  the  following 
formula. 


j  (E,Q)  = 


32  eV 


e  n  Ei  e  0 

7t  S^V2nm 


2e0 


(vi+Vgf  dEi 


The  quantity  vi  is  the  velocity  of  the  incident  ions  (ions  at  energy  Ej  incident  in 
the  look  direction).  The  quantity  Vg  is  the  velocity  of  the  detector  and  the  disk 
(negative  the  ram  velocity). 
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